Performance of a parasitic plant and its effects on hosts depends on the interactions between parasite seed family and host species

Abstract Root hemiparasitic plants act as keystone species influencing plant community composition through their differential suppression of host species. Their own performance also strongly depends on interactions with host species. However, little is known about the roles of parasite genetic variation vs. plasticity in these interactions. We grew plants from eight maternal families of the root hemiparasite Rhinanthus alectorolophus with six potential host species (two grasses, two legumes and two forbs) and without a host and measured fitness-related and morphological traits of the parasite, host biomass and overall productivity. Parasite biomass and other traits showed strong plastic variation in response to different host species, but were also affected by parasite maternal family. Parasite seed families responded differently to the hosts, indicating genetic variation that could serve as the basis for adaptation to different host plants. However, there were no negative correlations in the performance of families across different hosts, indicating that R. alectorolophus has plastic generalist genotypes and is not constrained in its use of different host species by trade-offs in performance. Parasite effects on host biomass (which may indicate virulence) and total productivity (host + parasite biomass) depended on the specific combination of parasite family and host species. Mean biomass of hosts with a parasite family and mean biomass of that family tended to be negatively correlated, suggesting selection for maximum resource extraction from the hosts. Specialization of generalist root hemiparasites may be restricted by a lack of trade-offs in performance across hosts, together with strong spatial and temporal variation in host species availability. The genetic variation in the effects on different hosts highlights the importance of genetic diversity of hemiparasites for their effects on plant community structure and productivity and for the success of using them to restore grassland diversity.


Introduction
Parasitic flowering plants represent c. 1 % of angiosperms (c. 4500 species, Heide-Jørgensen 2008). These species attach to the shoots or roots of their host plants via special organs (haustoria) and extract water, nutrients and carbon compounds from them (Těšitel et al. 2010). Many parasitic plants are hemiparasites which have green leaves and maintain the ability to photosynthesize . These species are generalist parasites which can grow with a wide range of species as hosts (Hautier et al. 2010;Matthies 2017), but their biomass (Matthies 2021), morphology (Campion-Bourget 1982;Jonstrup et al. 2016), patterns of allometry (Matthies 1995) and reproduction (Seel et al. 1993;Matthies 2017) can vary enormously depending on the host species. While some plant species are resistant against parasite infection and the growth of parasites with them is worse than without a host (Atsatt and Strong 1970;Cameron et al. 2006), others increase parasite performance strongly in comparison to autotrophic growth (Seel et al. 1993;Matthies 1995).
The interactions between parasitic plants and their hosts show many parallels with those between herbivorous insects and their host plants (Atsatt 1977;Pennings and Callaway 2002). However, while hemiparasites often have broad host ranges, many insects are specialist feeders on one or a few plant species (Jaenike 1990;Forister et al. 2015). A common explanation for this specialization of insects is the 'the jack of all trades is the master of none' principle, which states that specialists should outperform generalists on any specific host and that selection should therefore favour host specialization (Futuyma and Moreno 1988;Jaenike 1990;García-Robledo and Horvitz 2011). If there is a cost of adaptation to a specific host, a genotype adapted to one host is predicted to be less able to take advantage of another host species due to antagonistic pleiotropy. There should thus be trade-offs in the performance of herbivores with different hosts and a genotype that is well suited to one type of host plant is expected to have a poor performance on another host species, expressed as negative genetic correlations (Futuyma and Moreno 1988). However, empirical tests of such trade-offs in cross-host performance have found them only in some plant-herbivore systems (Henry et al. 2008;Gompert and Messina 2016), but not in others (Forister et al. 2007;Agosta and Klemens 2009).
While the genetics of the interactions between herbivorous insects and their host plants have been intensively studied (Forister et al. 2012), few studies have investigated the genetics of hemiparasite-host interactions (but see Mutikainen et al. 2000;Koskela et al. 2002;Ahonen et al. 2006). Individuals from different populations of hemiparasites grown in a common environment have been found to differ in a variety of traits (Zopfi 1993;Mutikainen et al. 2000;Rowntree et al. 2011). Some of this genetic variation in traits is likely to represent adaptation to different environmental conditions, e.g. habitat management (Zopfi 1993). Hemiparasite performance with individual host species has also been shown to differ among parasite populations and maternal families (Campion-Bourget 1982;Jonstrup et al. 2016;Sandner and Matthies 2017), indicating genetic variability in parasite responses to different host species, a precondition for the evolution of genotypes adapted to different hosts. However, the single study that has investigated possible trade-offs in the performance of hemiparasite genotypes across different hosts (Ahonen et al. 2006) did not find negative correlations between the biomass produced by the parasite Rhinanthus angustifolius with the grass Agrostis capillaris and the legume Trifolium pratense.
Root hemiparasites can have strong negative effects on their hosts (Ameloot et al. 2005;Phoenix and Press 2005;Press and Phoenix 2005;Matthies 2017). As host species vary in their sensitivity to hemiparasites, the parasites can change the competitive balance between plant species (Matthies 1996;Demey et al. 2015;Rowntree and Craig 2019) and act as ecosystem engineers. They often also reduce total productivity of the community (Davies et al. 1997;Mudrák and Lepš 2010;Borowicz et al. 2019). The outcome of the infection depends on the genetics of both parasites and hosts, as the impact of parasitism may differ between parasite populations (Mutikainen et al. 2000;Sandner and Matthies 2017) and seed families (Ahonen et al. 2006), and also depend on host genotype (Rowntree et al. 2011). However, it is not known whether hemiparasite seed families that benefit most from a host also cause the greatest damage to it, i.e. are most virulent.
To investigate the role of genetic variation in the interactions between a root hemiparasite and its hosts, we grew eight maternal families from a single population of the hemiparasite Rhinanthus alectorolophus with six potential host species (two grasses, two legumes and two forbs) in a common environment. Since R. alectorolophus is capable of autotrophic growth (Sandner and Matthies 2017), we also grew all parasite families without a host for comparison and to test for correlations between autotrophic and heterotrophic performance. A positive correlation between autotrophic and heterotrophic growth of seed families would suggest that some parasite seed families are generally more vigorous than others, while a negative relationship would suggest that an improved ability to extract resources from a host reduces the ability to grow autotrophically, increasing the degree of dependence on the host (Ahonen et al. 2006). The variation among families will indicate the evolutionary potential of the studied hemiparasite population, while the results of studies on the differentiation among hemiparasite populations (Campion-Bourget 1982;Mutikainen et al. 2000;Jonstrup et al. 2016;Sandner and Matthies 2017) reflected past genetic processes. Our study included a wider range of host species (six instead of two) than the related study by Ahonen et al. (2006), allowing us a more general test of the role of variation among families in the interactions between hemiparasites and hosts. The use of more hosts also allowed us to investigate whether those host families that are most beneficial for the parasites are also those that produce the least biomass with the parasites. However, due to restrictions on glasshouse space we had to use fewer parasitic families (8 instead of 25) than Ahonen et al. (2006).
We measured fitness-related traits and morphological characters of the parasite that have been used to delimit infraspecific taxa and ecotypes, host biomass and total productivity. We addressed the following specific questions: (i) Are traits of the parasite influenced by host species and maternal family, and do maternal families of R. alectorolophus react differently to individual host species? (ii) Are there trade-offs between the performance of parasite families with different hosts? (iii) Are the effects of R. alectorolophus on host biomass and total productivity influenced by parasite maternal family and are the negative effects of the different parasite families on the hosts and the performance of those families correlated?

Study species
Rhinanthus alectorolophus (Orobanchaceae) is an annual facultative hemiparasite growing up to 80 cm. It is usually pollinated by insects, but selfing is possible (Sandner and Matthies 2017). Rhinanthus alectorolophus is commonly found in Central European grasslands, especially in areas of low productivity and high light availability (Těšitel et al. 2015), and was in former times also a weed of cereal crops (Hartl 1974). It is one of the most frequently used species for studies on the ecology of hemiparasites, together with R. minor and R. angustifolius (Těšitel et al. 2011).
As host plants we selected two grasses (Dactylis glomerata and Lolium perenne), two legumes (Medicago sativa and Trifolium repens) and two non-leguminous forbs (Sanguisorba minor and Sinapis alba), which will be referred to by their genus names in the following. All these species occur together with R. alectorolophus.

The experiment
Seeds of R. alectorolophus were collected in a large population in northern Hessen (Germany) from eight large mother plants that were separated by at least 50 m. Seeds from the same mother plant (seed families) are at least half-sibs, but could also be full-sibs or be even more closely related because they resulted from self-fertilization. To break dormancy, the seeds were placed on moist filter paper in Petri dishes and kept at 5 °C until cotyledons had formed. Seeds of the host species were obtained from a commercial supplier (Appels Wilde Samen, Darmstadt, Germany) and germinated in Petri dishes at room temperature. Two seedlings of each host species were planted into 48 pots (11 cm × 11 cm) filled with a 4:1 mixture of commercial potting soil (TKS, Floragard, Interactions between seed families of a parasite and host species Oldenburg, Germany) and sand. After 2 weeks of growth, one seedling of the parasite was planted into each pot. For each combination of parasite maternal family and host species eight replicates were set up. In addition, five parasites from each family were planted into pots without a host. Plants were kept in a growth chamber at 16 h of light by metal halide lamps and a 20 °C/15 °C (day/night) temperature regime. Pots were watered regularly and their position randomized every 2 weeks. To follow the development of the parasites, the length of their longest leaf was measured after 2 and 4 weeks of growth and the day when a parasite produced its first flower was recorded.
After 9 weeks of growth, when the parasites were well developed and all but the smallest were flowering or fruiting, several traits were measured for each parasite: the height of each plant, the length of the vegetative part of the stem, total branch length as the sum of the length of all branches plus the height of a plant, the number of nodes below the lowermost flower (vegetative nodes), the length of the first five internodes, the number of flowers on the main inflorescence, the total number of flowers, the length of a fully developed flower, the diameter of a ripe fruit and the length and width of the longest leaf. Parasite leaf chlorophyll content was measured with a portable chlorophyll meter (SPAD-502, Minolta, Osaka, Japan) and the values obtained transformed into actual chlorophyll concentrations using the formula for total chlorophyll content given by Richardson et al. (2002). Parasites and hosts were harvested separately above ground, and all plant material dried for 48 h at 80°C and weighed.

Statistical analyses
Effects of host treatment and parasite family on parasite traits were studied by two-factor ANOVAs. We considered host treatment as a fixed factor and parasite maternal family as a random factor. According to the rules for mixed models (Zar 2010), the effect of host treatment was tested against the host by family interaction, while the other effects were tested against the residual variation. Parasite biomass, total branch length, the length of the first five internodes and the number of flowers were log-transformed to obtain normally distributed residuals. Simple main effects were calculated to test family effects within each host treatment. To investigate whether host species identity and maternal family influenced the allometric relationships between traits of the parasite and its biomass, we tested the effects of host and family adjusted for the effects of biomass in a general linear model that included the two factors and (log) parasite biomass as a covariate. The effects of biomass, host and the biomass by host interaction were tested against the interactions of these effects with family, while all other effects were tested against the residual.
To test if R. alectorolophus showed trade-offs in its performance with different host species or grown autotrophically, we related the mean biomass of the individuals of each parasite family when grown with a host to that produced with each of the other host species, and to that grown without a host. We studied the effects of host species and parasite maternal family on host biomass and total productivity (parasite + host biomass) by two-factor ANOVAs. To test for a possible relationship between parasite damage to a host and its performance, we related in a general linear model the mean biomass of each host grown with each parasite family to the mean biomass of that family grown with that host, after adjusting for the overall host effect. We presented this relationship in a partial residual plot created using the R-package visreg (Breheny and Burchett 2017). The relationship between the biomass of the individual parasites and that of their hosts was analysed by regressing the biomass of the host plant in each pot on that of the parasite separately for each host species. Statistical analyses were carried out with R 4.1.1 (R Core Team 2021).

Effects of host and maternal family on parasite biomass
The growth of R. alectorolophus as measured by the length of its longest leaf differed among the host treatments already after 2 weeks (F 6, 398 = 15.24, P < 0.001; see Supporting Information- Fig. S1) and these differences became more pronounced over time. The biomass of the parasite R. alectorolophus at the end of the experiment was also strongly influenced by the host treatment (Table 1). Lolium was the best host, followed by a group of species that were of a lesser quality as hosts (Trifolium, Medicago, Dactylis and Sanguisorba; Fig. 1A). The biomass of R. alectorolophus grown without a host was much lower, and that of parasites grown with Sinapis was even lower than without a host. The size of parasites grown with the best host Lolium was 21 times that of the size grown with the worst host Sinapis. The type of host treatment influenced parasite biomass even if the data for the very poor host Sinapis and for parasites growing without a host were omitted from the analysis (F 4, 28 = 10.7, P < 0.001).
The mean biomass of the parasites from the eight maternal families varied strongly and there was a significant interaction between the effects of host treatment and parasite family, i.e. the parasites from the individual families responded differently to the host treatments and the reaction norms of the different families crossed each other (Table 1; Fig. 1B). The analysis of simple main effects showed that parasite families grown with Lolium, Trifolium, Medicago Dactylis and Sinapis differed significantly in their biomass, while the variation among families grown with Sanguisorba or without a host was very small. Omitting the two worst host treatments did not change the results qualitatively (family: F 7, 268 = 11.1, P < 0.001; host × family: F 28, 268 = 1.7, P = 0.023).

Effects of host and maternal family on other parasite traits
Both host species and parasite maternal family significantly influenced not only parasite biomass, but also all other measured traits, and in most cases there was a significant interaction between the two factors (Table 1). Most traits were strongly correlated with biomass, i.e. parasite size (Table 2). Thus, similar patterns as for biomass were observed in most other traits in response to host treatment (Table 3) and maternal family [see Supporting Information- Fig. S2]. The strong correlation between the number of flowers and biomass indicated that biomass is a good proxy for fitness in R. alectorolophus. Trait values were usually highest for parasites grown with Lolium and lowest for those grown with Sinapis. Parasites with Lolium grew higher, needed less time to start flowering and produced more flowers and larger flowers and fruits than parasites grown with less beneficial hosts (Table 3). However, the correlations between parasite biomass and leaf chlorophyll content or length of the lowest five internodes were much weaker and there was only a very weak relation to the number of vegetative nodes (Table 2), which was more strongly influenced by maternal family than by host treatment (Table 1; Fig. 2A). Leaf chlorophyll content was highest for parasites grown with the legumes Medicago and Trifolium, and relatively low not only with Sinapis and without a host, but also with the good host Dactylis (Fig. 2B). Host treatment and maternal family significantly affected all studied traits of the parasite even after adjusting for differences in biomass, indicating that their effects on the various traits were not simple side effects of their effects on parasite growth, but that both host treatment and maternal family influenced the morphology and architecture of the parasite [see Supporting Information- Table S1].

Possible trade-offs between the performance of parasite families with different hosts
Negative correlations between the biomass of seed families grown with two host species would indicate trade-offs in the performance with different hosts. However, we did not find trade-offs between the biomass of families when growing with different hosts or without a host (   Hosts are in the order of decreasing mean parasite biomass. Significant differences among families within a host treatment: *P < 0.05; **P < 0.01; ***P < 0.001. Note log-scale for biomass. Interactions between seed families of a parasite and host species S4d). The positive relations were strongest between the biomass of the maternal parasite families grown with Dactylis and Medicago and Medicago and Trifolium (Fig. 3).

Host biomass and its relation to parasite biomass
The biomass of the hosts varied strongly among the six species (F 5, 35 = 17.04, P < 0.001), from 3011 mg for Lolium to 5540 mg for Dactylis, but depended also on parasite family (F 7, 320 = 3.42, P < 0.002) and its interaction with host treatment (F 35, 320 = 2.48, P < 0.001). Simple main effects showed that parasite family significantly affected the biomass of Dactylis, Lolium, Medicago and Sinapis, but not that of Sanguisorba and Trifolium (Fig. 4A). Total productivity (host + parasite biomass) also depended on host treatment (F 5, 35 = 9.90, P < 0.001) and the specific combination of host and parasite family (F 35, 320 = 2.56, P < 0.001), while the main effect of family was weak (F 7, 320 = 1.92, P = 0.067). The biomass of a host with a certain parasite family tended to be negatively associated with the biomass of that parasite family (P = 0.079; Fig. 4B). The biomass of individuals of the two grasses Dactylis and Lolium decreased significantly with the biomass of the parasite R. alectorolophus grown in the same pot ( Fig. 5A and B), which can be related to stronger effects of larger parasite individuals on their hosts or to larger hosts exerting stronger competitive effects on the parasite. In contrast, negative but non-significant relationships between host and parasite biomass were found for combinations of R. alectorolophus with Medicago, Sanguisorba or Trifolium (Fig.  5C, D and F). The biomass of Sinapis and R. alectorolophus (Fig. 5E) was positively related, but the parasites growing with Sinapis were very small.

Influence of the host plants on parasite performance
Five of the six host species (the grasses Lolium and Dactylis, the legumes Trifolium and Medicago, and the forb Sanguisorba) strongly increased growth and reproduction of the hemiparasite R. alectorolophus in comparison to that of parasites grown autotrophically, confirming that root hemiparasites like R. alectorolophus have a wide host range (Sandner and Matthies 2018;Matthies 2021). However, the performance of R. alectorolophus grown with the forb Sinapis was worse than that of parasites grown without a host, indicating resistance of Sinapis against parasitism or a low quality of the solutes the parasite extracted from this host. Negative effects of unsuitable host species on hemiparasites have been observed in other studies (Atsatt and Strong 1970;Ahonen et al. 2006;Cameron et al. 2006) and could be due to competition for light, water and nutrients by the host (Matthies 1995;Cameron et al. 2005), or related to the costs for the parasite of investing into non-functional haustoria (Ahonen et al. 2006). Grasses and legumes have been found to be generally good hosts for Rhinanthus spp. (Seel et al. 1993;Matthies 2021), but hosts from each functional group may provide hemiparasites with different compounds (Govier et al. 1967). Hemiparasites can obtain large amounts of carbon compounds from grasses (Těšitel et al. 2010), while legumes can be especially important for providing nitrogen due to their association with nitrogen-fixing bacteria (Govier et al. 1967;Haynes 2021). In our study, the grass Lolium was a significantly more beneficial host for R. alectorolophus than the other species. The two legumes, although not the best hosts in terms of biomass, resulted in a particularly high chlorophyll content of the parasites, presumably due to their provision of nitrogen. Leaf chlorophyll content of the related parasite Melampyrum arvense was strongly related to nitrogen content and higher in parasites grown with legumes than with hosts from other functional groups (Matthies 2017).

Genetic variation in parasite fitness with different host species
The performance of R. alectorolophus was strongly influenced by parasite maternal family, indicating genetic variation in growth and reproduction within the studied population, in contrast to the results of Hautier et al. (2010), but in line with among family variation found in R. angustifolius (Ahonen et al. 2006) and among population variation in other studies (Lammi et al. 1999;Mutikainen et al. 2000;Jonstrup et al. 2016;Sandner and Matthies 2017). However, plants from all the seed families could survive and reproduce successfully on several host plants, suggesting that they have plastic generalist genotypes (García-Robledo and Horvitz 2011).
Other factors apart from genetic variation (both additive and non-additive) could also have contributed to the observed differences among maternal seed families. The maternal environment can affect the quantity and quality of starch reserves, mRNAs, proteins, hormones and other metabolites packaged into seeds (Herman and Sultan 2011). However, these effects of the maternal environment on progeny are largest during the early stages of the development and diminish over time (Roach and Wulff 1987;Herman and Sultan 2011;Auge et al. 2017). Late traits like the biomass of adult R. alectorolophus are thus less likely to be affected by these maternal effects. Another potential mechanism for effects of the maternal environment on progeny are epigenetic changes, resulting in transgenerational plasticity (Galloway 2005;Herman and Sultan 2011). However, in experiments that found these epigenetic effects they were mostly induced by strong stresses (Hauser et al. 2011). In the current experiment, all seeds were sampled from large plants that had apparently not experienced very different stress levels.
The fitness of the parasites depended also on the specific combination of parasite family and host treatment (significant family by host treatment interaction). This variation in the performance of parasite families could serve as the basis for their adaptation to different host species. However, there were no negative correlations between the performance of the parasite families across different host treatment. Instead, we found in nearly all cases positive correlations, because the overall differences among parasite families in performance were strong. Negative genetic correlations would have indicated trade-offs  between the fitness of the parasites with different hosts due to antagonistic pleiotropy which are thought to strongly favour host specialization in parasites (Forister et al. 2007;García-Robledo and Horvitz 2011). The absence of tradeoffs indicates that adaptations increasing the fitness of R. alectorolophus with one host species are not likely to reduce its fitness when growing with other host species, reducing the likelihood of specialization to particular hosts. Our results agree with those of Ahonen et al. (2006), who found no tradeoffs but positive correlations between the performance of R. angustifolius across two host species, and several studies on herbivorous insects (e.g. Agosta and Klemens 2009; García-Robledo and Horvitz 2011; Laukkanen et al. 2013).
There are several possible reasons for the absence of tradeoffs in the performance of R. alectorolophus with different hosts and a lack of specialization. The typical habitats of R. alectorolophus, like that of most other root hemiparasites, are species-rich and thus for the parasites spatially heterogeneous (Těšitel et al. 2015;Holá et al. 2017), because the identity of host species available for an individual parasite is largely *P < 0.05; **P < 0.01; ***P < 0.001. The colour code for each family is the same as in Fig. 1. (B) Partial residual plot of the relationship between the mean biomass of each host species grown with each of the eight different seed families of the hemiparasite R. alectorolophus and the mean biomass produced by these families with that host, adjusted for the overall effect of individual host species (host species: F 5, 41 = 16.6, P < 0.001; parasite biomass: F 1, 41 = 3.2, P = 0.079). unpredictable. Such host environments that are variable in space and time favour the evolution of generalists (Futuyma and Moreno 1988;Jaenike 1990;Laukkanen et al. 2013). In contrast to many insect herbivores (Forister et al. 2015), transmission of root hemiparasites (dispersal) is not directed towards particular host species, which also creates variability in host species over time. If dispersal is not directed towards suitable hosts, specialization incurs high costs in parasites and frequent transmission among hosts is likely to lead to generalist parasites with suboptimal virulence (Rigaud et al. 2010). Root hemiparasite individuals can also attack several hosts simultaneously (Gibson and Watkinson 1989;Holá et al. 2017) and unspecialized individuals may thus benefit from a mixed diet (Marvier 1998;Sandner and Matthies 2018;but see Matthies 1996), because different host species may provide them with complementary resources (Govier et al. 1967;Sandner and Matthies 2018).
It could be interesting to study potential specialization on different hosts for parasite populations that grow in habitats in which single host species are strongly dominant. For example, parasite taxa that formerly grew as weeds of cereals like R. alectorolophus ssp. buccalis, R. angustifolius ssp. apterus and M. arvense had a far more predictable host environment that may have favoured specialization on a single host species. However, these hemiparasites growing as weeds of crops have strongly declined and have become threatened in many regions (Kornás 1988;Zając and Zając 2014).
There was no trade-off between autotrophic growth of the parasites and their performance with the six hosts, indicating that an increased capability of extracting resources from a host does not lead to a stronger dependency on a host (i.e. no tendency towards holoparasitism). An ability to grow without a host is important for root hemiparasites during the early stages of their life cycle when they have not yet access to a host and also in cases when the available hosts are unsuitable (Atsatt and Strong 1970).

Effects on other parasite traits
Both the host treatment and maternal family strongly influenced not only the biomass of the parasites, but all other traits. A large part of the variation in nearly all traits could be explained in terms of changes in parasite size, as most traits were strongly correlated with biomass. However, as in other hemiparasites (Matthies 2017), significant effects of host treatment remained even after adjusting for differences in biomass, i.e. parasite individuals of the same biomass but subjected to different host species varied in morphological traits like height, branch length, length of internodes, number and size of flowers, and fruit size. Host plant identity thus also influenced parasite allometry, architecture and morphology. Variation due to host species may partly be responsible for the extensive intraspecific morphological variation commonly observed in root hemiparasites (Zopfi 1993;Pleines et al. 2013) and suggests caution when trying to delimit infraspecific taxa on the basis of observational studies (Houston and Wolff 2012;Matthies 2017).
In contrast to other traits, the number of vegetative nodes was only weakly influenced by host treatment and much more strongly by maternal family, confirming that this trait is in hemiparasitic Orobanchaceae less plastic and under stronger genetic control (Campion-Bourget 1982;Zopfi 1993;Jonstrup et al. 2016). It is also a trait that varies consistently among seasonal ecotypes (Zopfi 1993). The time to flowering is negatively related to the number of vegetative nodes in Rhinanthus, but also influenced by host species (Jonstrup et al. 2016;Wesselingh 2016;Matthies 2021). In our experiment, flowering time of R. alectorolophus was mainly influenced by seed family, as variation due to the five good host species was very small. Our results support the idea that there is genetic variation within hemiparasitic populations that can allow them to quickly respond to changes in the environment, e.g. due to management, with changes in flowering time, resulting in seasonal ecotypes (Zopfi 1993).

Genetic variation in the effect of the parasite on the hosts
The biomass of the host plants varied more than 2-fold depending on parasite seed family, suggesting genetic variation in the effect of different families of R. alectorolophus on the hosts (i.e. their virulence). Genetic variation in parasite virulence has been found for R. angustifolius families grown with two hosts (Ahonen et al. 2006), among populations of R. angustifolius grown with A. capillaris (Mutikainen et al. 2000), and for R. alectorolophus grown with various host plants (Sandner and Matthies 2017). In contrast, there was no effect of population identity on the virulence of R. angustifolius or R. minor grown with Hordeum vulgare (Rowntree et al. 2011).
The effect of the parasite on total productivity also depended on the specific combination of host treatment and parasite family. This indicates that the success of using hemiparasites to reduce grassland productivity (Ameloot et al. 2005), promote plant diversity (Davies et al. 1997;Bullock and Pywell 2005), increase invertebrate abundance (Hartley et al. 2015) and control invasive plants (Těšitel et al. 2020) will depend on the genetic diversity of the parasite. Previous research already found that genetic diversity appeared to promote establishment of the related parasite R. minor in a high species environment (Rowntree and Craig 2019).
In this study, parasite families that across all host species produced the highest biomass tended to be associated with the smallest hosts, suggesting that they have the strongest negative effects on the hosts. There may thus be selection for increased negative effects on hosts. This is in contrast to the optimum degree of virulence and the prudent rather than maximum resource extraction predicted by the model of Hautier et al. (2010) for hemiparasites. However, these authors also cautioned that if several parasite individuals parasitize the same host individual, maximum resource extraction would be favoured (see also Regoes et al. 2000;Rigaud et al. 2010). Maximum use of host resources by hemiparasites that strongly reduces host growth might also be selected because the hosts are not only an important source of nutrients, carbon and water for the hemiparasites (Těšitel et al. 2011), but also competitors for light (Matthies 1995).
The biomass of the parasite and that of the grasses (both Lolium and Dactylis) grown in the same pot were significantly negatively correlated. This could indicate that large parasites suppressed hosts more strongly than small ones, but also that large host individuals were stronger competitors and prevented strong growth of the parasites. However, similarly strong negative correlations between parasite and host biomass were not found for the other host treatments indicating no clear relationships between resource extraction and parasite vigour. In the case of Medicago and Trifolium this may have been due to a tolerance of parasitism, as these species are good hosts for R. alectorolophus, but little damaged by parasitism (Matthies 2021).
Surprisingly, there was a significant positive relationship between the size of parasites and Sinapis individuals grown in the same pot. As Sinapis was clearly an unsuitable host one would have expected that the negative effects of the Sinapis plants on R. alectorolophus due to competition for light and nutrients would increase with their size. A possible explanation for the positive correlation is that the parasites may nevertheless have obtained some resources by parasitism and this benefit may have been positively related to individual host size.

Conclusions
We found significant variation among seed families from a population of the hemiparasite R. alectorolophus in their performance with different hosts. This indicates evolutionary potential and could provide the basis for specialization to different host species. However, absence of trade-offs between the performance with different hosts support the notion that hemiparasites, like many herbivorous insects (García-Robledo and Horvitz 2011), are not constrained in their use of different host species by trade-offs in performance. Together with the strong spatial and temporal variation in host species availability and gene flow, the lack of trade-offs could explain the absence of host specialization within the parasite population (Rigaud et al. 2010). Increased damage to hosts (virulence) tended to be positively correlated with parasite fitness in R. alectorolophus, suggesting selection for maximum resource extraction from hosts. The variation among seed families of R. alectorolophus in their effect on different hosts highlights the importance of the genetic diversity of hemiparasites for their effects on the community structure and diversity of grasslands.

Supporting Information
The following additional information is available in the online version of this article. Table S1. Effect of host species (H), maternal family (F), parasite biomass (B) and their interactions on traits of the hemiparasite R. alectorolophus. p < 0.05 are in bold face. Table S2. Trait means for the different parasite maternal families. Different letters in the exponent denote differences between families (Tukey-test, P < 0.05). Figure S1. Growth of R. alectorolophus (length of longest leaf) with six different host species and without a host (Lol, Lolium perenne; Tri, Trifolium repens; Med, Medicago sativa; Dac, Dactylis glomerata; San, Sanguisorba minor; no, no host; Sin, Sinapis alba). Vertical lines show ± 1 SE. Figure S2. Reaction norms of several traits of eight maternal families of the parasite Rhinanthus alectorolophus in response to the host treatments (six different hosts and no host). (a) Height, (b) length of the vegetative part of the stem, (c) total branch length, (d) length of the first five internodes, (e) days until flowering, (f) number of flowers on the main inflorescence, (g) total number of flowers, (h) length of flowers, (i) diameter of fruits, (j) length of the longest leaf, and (k) width of the longest leaf. The colour of each family is the same as in Fig. 1. Host species are in order of decreasing parasite biomass. Significant differences among families within each host: *, p < 0.05; **, p < 0.01; ***, p < 0.001. For abbreviations of host species names see Fig. S1. Figure S3. Pairwise relationships between the mean biomass of individuals of the hemiparasite R. alectorolophus belonging to eight different families grown with different hosts. Continuous regression lines are shown if p < 0.05, dashed lines if 0.05 < p < 0.1. Figure S4. The relationship between the mean biomass of individuals of the hemiparasite R. alectorolophus belonging to eight different families when grown with different hosts and without a host. A dashed regression line is shown if p < 0.1.